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ABSTRACT 


Recently , a fluid turbine which has a ^ '’SS 

designed and tested with wa . a tube ^ The Montreal 

turbine-brush is to clean up Products from refrigeration 

Protocol prohibits the use of CFC products^ * n 1996 . 

industry or from ^ ndus 1 t i 1 ^ ers 9 devices or a combination of 
Alternatives for the ^ 1 ® ans ® ' r sho uld be found. One of the 

alternative devices with a cleanser snouiao fouling with a 

methods is to develop a evice describe a turbine connected 

cleaning medium. In this P *P rbine with the brush should be simple 
with a brush. However , the turbine wi Lth ‘ the small liquid 

and easy to install^ This d is activated by the liquid 

turbine with a brush, ine h - turns the brush cleaning 

flowing through the tube. Then the t^rb r n and the 

fouling along the tube. Based ^fricai relationship of drag 
Bernoulli equation along W1 th ^ D Daratus, a relationship of 

force obtained from an experi and aeometric variables 

the rotational speed, the nmber of ™?Sonal speeds 

£ -^^^^nMSons. Further worK was 
recommended for improvements . 


SUMMARY 


One of methods of cleaning fouling in a tube will be a turbine- 
brush unit to replace R113 as a cleanser after 1995. The turbine- 
brush consists of a turbine and a brush with a connector between 
the two. The turbine-brush was created in KSC,NASA and demonstrated 
its rotational motion while it was held at a position. 

Objective of the program is to optimize the geometry of the 
turbine-brush including the blade angle, the number of blades, and 
arrangements of the connector and brush, flow rates, rotational 
speeds, diameters of tubing and performance of cleaning before a 
final product is designed and manufactured. 

As an initial step, a relationship among the above variables was 
developed based on the energy conservation law and the force 
momentum theorem. Numerical values of an existing turbine-brush are 
used with the relationship to predict the rotational speed. 
Predicted rotational speeds are compared with observed rotational 
speeds. These two differ in one order of magnitude. 

The drag coefficients were computed based on the incoming velocity . 
From the open literature, we did not find any information of its 
drag force coefficients, because this may be the first device to do 
the cleaning in the tube with the fluid as driving energy source. 
Possible alternative coefficients useable for the turbine-brush are 
those for a spinning sphere or a disk and cylinder connected in 
tandem. From the experimental results, the force was measured and 
its coefficient was calculated. As a result, the equation relating 
the variables produced the angular velocity. 

While it rotates, the flow downstream of the turbine-brush is a 
combination of a usual turbulent one dimensional flow with a forced 
vortex motion due to the turbine motion within the connector and 
brush then one more complicated form of forced vortex motion takes 
place downstream of the brush. A laminar sublayer flow was assumed 
to pass through the gap between the turbine-brush and the wall of 
the tubing. 

A simple numerical model was developed. One hundred sixty eight 
volume elements were constructed in order to use in a three 
dimensional unsteady flow conditions. Boundary conditions, 
properties , turbulence model, and so on have to be supplied to run 
and yield a result. 

Further work is recommended as follows: (l)a coded program of 
computation is recommended for various parameters of the turbine- 
brush, ( 2 )measurement of drag force coefficients in the wide range 
of the volume rates, (3 ) study of constraints under which the 
relationship derived be optimized including the degree of the 
reaction of the turbine, ( 4 ) experiments be performed with a 
different forms of the brush,(5)a numerical simulation be continued 
with the model to determine a feasibility of the numerical scheme 
for the turbine-brush in the future. 
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NOMENCLATURE 


A 

B 

C D 

d 

F 

g c 

H 

h 

K.E. 

m f 

P 

P . E . 

Q 

r 

u 

U 

V 

v 

W 


Area of tube cross-section, in , 

Width of the brush, or the length of turbine-brush 

sections , m 
Drag coefficient 
Diameter of the tube, m 


Fores N ■ 

Constant defined by the force: unity is the SI unit 

Enthalpy, kJ/s 

specific enthalpy, kJ/kg 

Kinetic energy, kJ/kg 

Mass rate, kg/s 

Static pressure, N/m z 

Potential energy, kJ/kg 

Heat transfer rate, J/s 

radii of the rotating elements, m 

Specific internal energy, kJ/kg 

Tangential Velocity, m/s 

Fluid velocity, m/s 

Specific volume, m 3 /kg 

Rate of work, kJ/s 


Greek Letters 


a 

p 

31 

P 

2 

CO 

p 


Turbine angle 

Weight adjusting factor 

3.14 

Density, kg/m 3 
Summation 

Angular velocity, /s 
Dynamic viscosity, N s/m 


Subscripts 


1-2 

1-4 

b.b. 

b.h. 

br 

D 

f 

in 

out 

shaft 

skin. f . 

rot . 

tension 

tur 


representing the length between p L and p 2 

representing the length between p x and p 2 

Brush Base 

brush and its hub 

Brush 

Drag force 

Fluid 

Entering section 
exiting section 
shaft work 
Shear force 

Referring to rotational parameters 

Tension of a wire 

turbine 


I . INTRODUCTION 


CFC refrigerants/cleansers will not be allowed to use in cleaning 
fouling in tubes starting in 1996. In this paper, we describe a 
device which will clean fouling in tubes. The device consists of 
two parts; a turbine and a brush with a connecting part. Details 
of the turbine, brush and the connector are shown in 
Fig. l(Werlink, 1993 ) . The turbine has four blades with 70 degree 
wide and there is a spacing of 20 degrees between the blades. These 
blades are installed in a direction of 45 degree from the axial 
direction. A connector links the turbine with a brush. 

The brush is made of nylon fibers so very densely structured that 
fluid may not pass through. It is possible to consider that the 
friction in the gap between the turbine and the tube wall and brush 
and wall may be of a nature of a laminar flow. However, the drag 
force of the turbine-brush unit is very complicated. At a glance, 
the unit may be assumed to be a spinning sphere (Goldstein) or a 
disk connected to a cylinder in tandem( Blevins ) . While a cleanser 
water solution passes through a tube, more or less, the friction 
force may be reduced due to the fact that the cleanser is a 
material made of long chain structured polymers (Sellin and Moses). 

On the other hand, a constant vortex motion provided by the turbine 
followed by a swirling motion of the fluid within the connector and 
the brush, again, a vortex motion due to the brush rotation are very 
complicated and difficult to analyze exactly. 

The combined vortex and turbulent flow may be analyzed by 
superimposing two separate flow characteristics together, i.e. , an 
axial turbulent flow passing through the geometry of the turbine- 
brush and a rotational flow within the turbine-brush configuration 
may be analyzed separately and results should be superimposed. 

In this project, preliminary results were obtained based on 
simplifying assumptions and measurements of drag forces. Based on 
the assumptions, the energy conservation law was used along with 
the force momentum theorem tieing up with the drag force in terms - 
of pressure drops across the turbine— brush. By doing so, the speed 
of the unit in rotational motion was predicted and compared with 
the experimental results . 


I I. ANALYSIS 


2.1 Drag Force and its Coefficient 

Consider a control volume shown in Fig. 2. Certain assumptions are 
needed to deal with the flow phenomena for a first order 
approximation solution. The assumptions used for the fluid 
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are (Munson et al.) s 


1 . One-dimensional incompressible steady flow 

2 . no heat transfer from or to the flow, and t lB W* 

3 Fully developed turbulent flow 

4 * The tube has a uniform internal diameter. . 

in computing the drag coefficient. 


Based on the force-momentum theorem, 
horizontal direction is the same 
change ( Munson , et . al ) : 


the sum of forces in the 
as the rate of momentum 


2 F =( m f / g c ) ( V out ) 

The sum of forces are 
forces, the drag force 


( 1 ) 

collection of the tension force, static 
and the skin friction forces (Fig. 2) : 


- F. 


tension 


+ (p L -p 2 )A + F D - F skin ,.i_ 2 . = ( m f/9c)(^2 " V l )(2) 


the presrrs ?"he 

stations 1 and 4 are used equation(2) becomes (Fig. ). 


— F tension (Pi 


p 4 )A + F D + F Bkin .f.i_4. 


(m f /g c )(V 4 - VJ = 0 


because V 4 is the same 
station 4 in the tube. 


as V L for a fully developed flow at the 
From the above, we obtain F 0 : 


F n = F 


tension 


+ (P 4 “Pl) A “ ^"skin.f .1-4. 


By the definition of the drag coefficient, C D , the drag force is 
given by: 


F d = p f V 1 2 AC D /(2g c ) 


( 3 ) 


C D = [ F tenBioI i +(P4 - Pl) A " F Bkin.£.l-4.]/(Pf V l 2A / 2 gc) 

The coefficient, C„ is computed from ^rictL^or^is 

tubing w£^!l .^Th^f rictioi^'f^>rc^ i^the^xia ] 3 direction has a form, 


Fakin.f .1-4. 




( 4 ) 
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2.2 Energy Conservation Law 


Consider the energy conservation on the control volume shown in 
Fig. 2 (Moran and Shapiro): 

Q - W = (K.E. + P.E. + H ) out - (K.E. + P.E. + H) ln (5) 

where Q is the heat transfer rate, W, the mechanical work and K.E. 
and P.E. are the kinetic and potential energy , respectively and H is 
enthalpy defined by: 


H= m f h ; h = u + p/p, m £ = mass rate 


( 6 ) 


Applying the usual assumptions to liquid steady flows,i.e., no heat 
transfer, temperature being unchanged, i.e, t ln = t out , and the 
density being constant, the above equation is written as follows: 

-W = (m f /g c ) [ (V 2 2 /2 + p 2 v) - (V 1 z /2 - p,v) ] (7) 

The rate of work in equation(7) includes the shaft work and the 
friction work as the unit rotates: 


W ^shaft ^skin.f. 


(8) 


Since the rotational force comes from the motion of the blades , 
all of energy of fluid is assumed to turn the turbine-brush, the 
shaft work will have a form: 

W 8haft = Jtot« 2 /2 ( 9 ) 

where co is the angular velocity and J tot is the sum of the moment of 
inertia about the axial direction. The sum of the moments consists 
of moments of inertia of the turbine, J tur , the connectors , J conn . f 
brush hub, J b . h . and the brush, J b: 


J tur * (pt) tur (n/2) (r 4 4 - r x 4 ) (cx/360)N|3 tur 

*iconn. = ( pi ) conn, (^t/2) (r 3 4 - T 2 ) Pconn. 

Jfc.h. = ( P® ) b.h. ( ) ( r 3 - r 2 ) b.h . 

Jb = (pB) b (7t/2)(r 3 4 - r 2 4 ) b 
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where p is the density, t is the thickness, a is the angle of the 
blade, and p is a number adjusting the geometry with respect to 
one of the components of the turbine-brush, N,the number of blades, 
and r's are radii, respectively (See details in derivations of 
APPENDIX) . 

Since the turbine-brush unit is held in a position but rotates, the 
friction work due to the translation motion is negligible. The 
friction between the brush and surface of the tube wall is due to 
the force normal to the centrifugal force. This force should be 
multiplied by the tangential velocity of rotation (See details in 
derivation of APPENDIX) : 


w, 


br , rot 


= (2n/3gJp dry (pB) br co 3 r 7 (r 7 3 - r 2 3 ) bu P b . h . 


(14) 


The magnitude of the dry friction coefficient between the turbine 
the brush and the tube wall is 0 . 35 (Baumeister and others ). Since 
the flow in the gap between the tube wall and turbine housing 
laminar, the work against the rotational motion is: 


W rot . = ( 2jt/3 g c ) fx f ( pB ) tur co 3 Nr 4 ( r 4 3 - r, 3 ) 


turWtur 


(15) 


After collecting all of these terms, equations ( 2 ),( 3 ),( 4 ),( 9-15 ) 
and substituting into equation(7), a relationship .?t 
variables such as the angular velocity, the speed of the unit along 
the tube, and the number of the blades is shown as follows. 


c 1 r 7 (r 7 3 - r 3 3 )o> 3 + c 2 r 4 (r 4 3 -r^JNto 3 - c 3 (r 4 4 - r^JNco 2 
( c 4 + c 5 ) (r 3 4 - r 2 4 )w 2 -c 6 (r 7 4 - r 3 4 )co 2 

+ c 7 (r 6 + r t )W -c 8 = 0 

where 

Cl = (2n/3g c )p dry (pB) br . Pb.h. 

C 2 = ( 2ll/3g c ) )lf ( pB ) turPtur ( o/360 ) 

c 3 = (Jt/2)(pB) tur P tur (a/360)/2g c 

c 4 = (ji/ 2) (pt) conn. Pconn. /2g c 

c 5 = (3t/2) (pB) b . h ./2g c 
c 6 = (ji/2) ( pB) b /2g c 
c 7 = m f /(8*g c ) 

c 8 = — ( m f / 2g c ) V 1 Z C D tFtgngionVj 


(16) 

(17a) 
(17b) 
(17c) 
(17d) 
(17e) 
(17f) 
(17g) 
( 17h) 
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where p b . h . is a factor adjusting the geometry for the brush and hub, 
Pconn t a factor adjusting the geometry for the connector, and p tur ^is 
a factor adjusting the geometry for the turbine. ^2 ^2tran8. 
V 2rot . 2 , v, 2 « v 2trans . 2 V 2rot = ( ( r 6 + ri )/2)0> were also used in 

equation ( 16) . Subscripts 1, and 2, here refer to the inlet and 
outlet , respectively . 


2.3 Application of the Derived Equation 


For the number of blades, N= 4, r's being read from Fig. 1, listed in 
Table 1, equation (16) becomes a third degree polynomial in to at a 
given volume rate,.i.e., V x = constant: 

a 3 (o 3 + a 2 co 2 - a 0 = 0 

where a's are expressed as follows: 
a 3 = c 1 r 7 (r 7 3 - r 3 3 ) + c 2 r 4 (r 4 3 - r x 3 )N 

a 2 = c 7 (r 6 + r x ) 2 - [c 3 N(r 4 4 - r x 4 ) + (c 4 + c 5 )(r 3 4 - r 2 4 ) 

+ c € (r 7 4 - r 3 4 ) ] 

a 0 = ^8 

where C D is assumed to depend on the incoming velocity, V x and 
constant determined from the experiments. Table 1 summarizes all of 
the numbers in the above computations. 


(18) 

(19a) 

(19b) 

(19c) 


Table 1 Values of material properties , geometry 
and flow properties 


Turbine 
r 3 =0. 00381m 


r 4 =0. 01905 
r 5 =0. 02159 m 
r 6 =0. 02301 m 
r 7 =0.0254 m 
t=0. 0015875 m 

B=0. 0047625 m 

p=7900 kg/m 3 


Brush 

r 2 =0 .00508 
r 3 =0 . 00762 

r 7 =0.0254 m 

B=0. 008382 m 
p=1121 . 625 kg/m 3 


Brush Hub 

r 2 =0 . 00508 
r 3 =0. 00762 

B=0. 008382 m 
p=4800 kg/m 3 


Connector 

r 2 =0. 00508 
r 3 =0 . 00762 

t=0. 0015875 
B=0. 00254 
p=7900 kg/m 3 
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The density and viscosity of fluid are p=997 kg/m 3 and p=885xl0 
Ns/m 2 . 


2.4 Numerical Simulation 


a software NEKTON ( Fluent , Inc . ) was made available on a DEC 

workstation “the branch through the KSC % 

a finite element based fluid problem solver that is capable of 
analyzing flows of three dimensional unsteady nature. This is 
result of a spectral element method developed by Professor Pa era 
of MIT and commercialized by Fluent, Inc. 

The flow domain of the turbine-brush is a cylinder with a length of 
Ud d being the diameter of the cylinder. At the location of 3d 
from the entrance, the turbine-brush of 2d long is placed. T 
location of the exit is then 6 d away downstream 

brush which may be the minimum distance for which a fully develop 
flow at the exit is warranted. 

Startinq from the entrance, 21 surface elements were constructed. 
These elements were moved up to the next position a 
cylinder length, thereby, became volume elements. At . the *°^ent 

of the turbine brush, the surface was . . dlv ^* d ni . that these 
concentric elements from the previous fluid elements so that these 
elements are identified as rotating elements at given locations. 
The connector and brush sections were made into volume elements 
a similar manner to the turbine section (Fig. 3 ) . 

At present, 168 finite volume elements were constructed for the 
flow domain and stored in the SUN workstation network at UNC- 
charlotte. Additional variables to be supplied are houn y 

conditions, fluid properties, rotating boundary ea ch "element' 
turbulence model, and deformation information for each elemen . 
This is ^ very crude model for the turbine-brush. However, this 
would yield important an insight so that an analyzer can determine 
whether a numerical study may be of help m the turbine brush 
development . 


III. EXPERIMENTAL INVESTIGATION 


3.1 Type of Probes and Their Locations 


Fig. 4 shows an experimental apparatus. A hydrant with ^ control 
valve suDDlies water to the test section of a transparent tubing of 
I d 2 in! Thl valve controls the flow rate. The turbine-brush unit 
is held bv a wire and a load cell transducer measures the tension 
force of the wire. The pressure ports are shown as p,, Eh. Pi' and p.. 
The rotation speed w!s measured by picking up pulses from a 
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rotating magnet which is placed in a turbine blade. The sampling 
rate was 49.98 per second. The load cell and pressure transducers 
have the measuring ranges as shown in Table 2. 


Table 2 Transducer Ranges 


Transducers 
Load cell 


Pi 

P 2 

Pi 
P 4 


ranges 

0-50 lb 

14.7-60 psia 
14.7 - 50psia 
12.0 - 30psia 
12-30 psia 


Accuracy in F.S. 

0.1 % 

0 . 02 % 

0 . 02 % 

0 . 02 % 

0 . 02 % 


A PC and labview software display and record the data in an ASCII 
file in a monitor and disk. A typical data set is shown in Table 3. 


3.2 Data from Experiments and V If F D Computation 


Applying the force-momentum theorem to the turbine-brush held at a 
point produces an equation in the form: 


“Ftenelon + ( Pi “ P 2 ) A + F Drag “ F skin.f.l- 2 . 

= (m £ /g c )(V 2 - V x ) (2) 


where the first two quantities of the left hand side of the 
equality and V L are known and the drag force is given by F Dra „ - 
C^OfV, 2 A/2q ( ., A being the cross section area of the tube. From this 
equation V 2 is obtained. Computation yielded that V 2 's are within 
a range of 95 percent of V x . Table 3 summarizes data from the 
experiments, computed velocities and the drag forces. The 
quantities in parentheses are the numbers in the Pound-Inch system. 

The reading devices have ranges of uncertainties within 0.10 
percent, yielding an uncertaintity of about 1.0 percent for the 
measurement system. 
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Table 3 Summary 


of Flow Rates, Pressures, Tension Force 
and Drag Force 


Flow Rate p L 

mV s 

kpa 

(gpm) 

(psi) 

.0037 

21.86 

(58.5) 

(3.17) 

.0047 

21.79 

(74.47) 

(3.16) 

.0052 

20.80 

(82.00) 

(3.02) 


Pz P 3 


kpa 

kpa 

(psi) 

(psi) 

10.34 

1.014 

(1.50) 

(0.147 

10.76 

1.048 

(1.56) 

(0.152 

10.74 

1.50 

(1.60) 

(0.218 


P< ^"tension 


kpa 

(psi) 

N 

(lbs) 

2.64 

39.44 

(0.38) 

(8.86) 

2.64 

41.15 

(0.38) 

(9.25) 

2.81 

39.41 

(0.41) 

(8.86) 


Drag Force 


m/ s 

N 

(ft/s) 

(lbs) 

1.822 

0.468 

(5.98) 

(.105) 

2.32 

0.53 

(7.61) 

(2.36) 

2.01 

0.655 

(6.58) 

(2.97) 


IV. RESULTS AND DISCUSSION 


4.1 The Angular Velocity 

in the computation procedure, the geometry factors were used; fW 
uom putarr ^ _ 0 . 65 . Table 4 shows the computed and 

observed f reciuencies for three cases * 

The predicted and the observed frequencies do not agree well. Jh 
predicted* 3 values are of magnitude of one order high we ' thar the 
observed values. This means that the mass of the rotating y 
has Sin estimated less than it should be and that the sum of the 
SrS of IXtfng parts of the turbine-brush should be estimated 
more carefully, rather than using the geometry factors. 

n _ p 0 f t he assumptions used was the flow in the gap between the 
turbine-brush and the tube wall is a laminar sublayer. This may no 
be true Estimating the frictional work in the rotating motion of 
the turbine-brush may not be accurate, particularly with i the _brush. 
More data are needed for the drag coefficients in a wide range of 
thf volume rates to know whether the deviation between the two 
va?uJs of thf angular velocities would be greater m the high 
Reynolds number side . 


4 . 2. Coefficients of Drag Forces 


326 


From Table 3, the drag coefficients were computed for three 
different flow rates. The Reynolds numbers and the coefficients are 
listed in Table 4. Further computed results will be compiled for a 
later report. 


Table 4 Coefficients of Drag forces at Flow Rates 


Re 

c D 

f pred. 

f observed 

1 . 08xl0 5 

0.14 

43.98 

7.0 

1.37x10 s 

0.43 

47.41 

7.0 

1.19x10 s 

0.20 

47.86 

7.0 


Flow Rate 
m 3 /s (gpm) 

0.0037(58.5) 

0.0047(74.47) 

0.0052(82.00) 

The difference between the incoming velocity and the moving speed 
of the turbine-brush is very little. The drag coefficient chang 
from 0.14 to 0.4 in the Reynolds number range at iO. The 
coefficient of the turbine-brush is compared to those of the 
spinning sphere (Goldstein) , and the disk connected to a cy in e 
?n ?anSL. The drag coefficients of the turbine-brush rs less than 
the spinning sphere and the disk connected to a cylinder (Blevins ) . 

4.3 Constraints under which optimum values are obtained 

An important constraint among others seems to be the degree of the 
reaction of the turbine-brush, which is defined as ( Shepherd ) : 


R = [(U, 2 - U 2 2 ) + (V r2 2 - V rl 2 ) ]/(2g c E) 


( 20 ) 


where E is defined as : 

E = [(V 2 - V 2 2 ) + (U x 2 - U 2 2 ) + (V r2 2 - V rl 2 ) ] / ( 2g c ) 


( 21 ) 


where all of the symbols used are denoted in Fig. 5. The absolute 
velocity V is resolved into a component V m passing through the 
axis ( I radial velocity) and a tangential component V . The 
reaction is an indication of energy transfer from the fluid to th ® 
turbine. For example, R = 0 implies that there is no jellied 
static head or pressure in the rotor and such a mach ^® 1S . 
an impulse type. In general, the term reaction alone is used 
whenever the turbine is not purely impulse, but for the |P e ^ a f 
case of steam turbines, it has come to imply 50 percent o 
degree of the reaction. The turbine of the turbine-brush has a 
deqree of the reaction between 0 and 50 percent. This aspect of th 
degree of reaction for the turbine-brush should be studied further. 
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Other important constraints are the optimum radii of ^structure 
W -, KP thpse determine the weight _ of the turbine Drusn^, 
consequently S the rotating inertia. The structure of the connector 

le^ms to be y a very important constraint. A " e " AheTTariab^s 
structure of the brush will be important. All of these variables 
are related with the geometry factors. The more accurate these are, 
the better predictions seem to be made. 

Strictly speaking, the actual application of this device in 
cluing fouling, the device will travel along the tube, not held 
at a place. The current experimental apparatus holds the turbine- 
brush at a point, thus the drag coefficient might be different from 
those obtained from the apparatus in which the turbine-brush was 
held at a point. A series of this type of experiments shoul 
carried out in the future. 


4.4 Numerical Simulation 


The numerical analysis be continued to establish a feasibili y 
study. This type of the feasibility study may be needed to justify 
the required time for a coding a new program or justify using 
commercially available code, if needed. 


V . RECOMMENDATION 


Further work is recommended as follows. 

1 Determine the weights of the turbine-brush as a unit and 
individual components. This will help determine the > mas .. . 
inertia of the turbine-brush and mass adjusting factors, because 
these control the prediction of the rotation speed. 

2 More drag force measurements may be needed from low volume 

rates. Additional experiments would tell us Wh ®f h ?L flows ' 5 
is reliable in the low or high velocity side of the flows. 

3 Investigation of the constraints under which the relationship 
derivedbe optimized including the degree "ft''®' 1 ”"”' 1 
turbine. Beciuse of the brush attached to the turbine, this may be 

a brand new study. 

4 Different forms of brushes be attached and experimented : for a 
^ssib?freductTon of the drag force. Rather than the radial brush, 
attach fibers in the connector along the axial direction s° 

, h fibers aet stretched due to the centrifugal force in rotation. 
This would leduce the drag force and enhance the rotational speeds. 
This kind of fiber arrangement could be beneficial for a smalle 

diameter . 
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5. Computation procedures be coded for the current verification of 
the equation as well as its future use. 

6. Numerical analysis with constructed elements be carried out to 
determine feasibility for further numerical study and to know 
whether it could aid the design development activity based on the 
software NEKTON or a different code or coding a new code. 


VI. CONCLUDING REMARKS 


Based on the Newton's law of force and the energy conservation law, 
an equation was derived relating certain number of geometrical 
variables with flow variables together. Pressure drops across the 
turbine-brush were measured in order to compute the drag force 
coefficients. These coefficients were used in the derived equation. 
This equation provided the rotations per second. These have been 
compared with the experimental values. They do not agree at all. 
The reason was discussed in the Results and Discussion section. 
Further work has been suggested in the Recommendation section. 


VI I. APPENDIX 

A. DERIVATION OF THE KINETIC ENERGY OF THE BLADES 

From Figs. 6 (a) and 7 (a), the mass moment of inertia of the arc 
element is : 

dl = r 2 (dm); dm=p(2nrdr)t=2ji(pt)rdr; dl=2jt(pt)r 3 dr 

I = J2ji(pt/g c )r 3 dr = (ji/2g c ) (pt) (r 4 4 - r^) (Al) 

where p is the density, and t is the thickness of a plate. 

Integration of the above with respect to the radius for a blade 
that has a degree, the outer and inner radii as r 2 and r 4 gives: 

I tur = (Jt/2g c ) (pt) tur (r 2 4 - r, 4 ) (a/360) (A2) 

Thus, the kinetic energy of the N turbine blade is: 

(jt/2g c ) ( pt ) tur co 2 ( r 2 4 - r^N/2 (A3) 
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B 


.DERIVATION OF THE OF THE ROTATIONAL FRICTIONAL FORCE AND ITS WORK 


From 
is : 


Figs . 6 (a) and 7(b), the centrifugal force of the arc element 


dF cent = rco 2 (dm); dm=p(2jirdr)B; w = constant 


(Bl) 


where u> is the angular velocity of the blade. In ^ration of the 
force with respect to the blade angle and radius with the number o 


blades gives: 

F C ent = ( 2n/3g c ) ( pt ) tur w 2 ( r 4 3 - r 1 3 )(a/360)N 


(B2) 


The friction force is, then: 

F frict = l^dry ^ 1 cent, 

where is the dry friction coefficient between the brush and the 
tublwali Thus, the rate of the frictional work is: 


W br ,rot = ( 2 Jt / 3 g c ) n dry ( pB ) br i*> 3 r 2 ( r 7 3 - r 2 ) 


(B3) 


C.THE NEWTON-RAPHSON METHOD FOR FINDING ROOTS OF A POLYNOMIAL 

Consider a function of x, f(x)=a 0 + a 2 x + a 2 x 

+a 3 x 3 + a 4 x 4 +• • ■ (Cl ) 

Suppose an estimated root is x„. Then evaluate the function 
f (x 0 ) , that is: 


f (x 0 ) = a 0 + a 3 x 0 + a 2 x 0 2 + a 3 x 0 3 + a 4 x 0 4 + 


(C2) 


Take the derivative of the function with respective to x and 
substitute x with x 0 : 

df(x 0 )/dx 0 = a, + 2a 2 x 0 + 3a 3 x 0 2 + 4a 4 x 0 3 + (C3> 

Now make a ratio of f (x 0 ) /df (x 0 ) /dx, subtract from xO and set it 
as : 


x n = x 0 - f(x 0 )/df(x 0 )/dx 


(C4) 
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Let us compare x n with x 0 in terms of their magnitudes ; hence the 
absolute value of the difference is less than a certain number , e 
=0.01 for which the accuracy is warranted: 


|x n - x 0 | <; e 


( C5 ) 


then, one of the roots of equation(8) is x n . For other roots, the 
above procedure should be repeated for the polynomial ( Stoecker ) . 

Using velocities of Table 3 and C D 's of Table 4 for three cases of 
experimental results, the roots were found. w= 276.2072, 297.7274, 
and 300.5572, respectively. The other two roots are imaginary 
numbers for the three cases. Thus co = 2jc(rotation/s ) produces 
43.98,47.41, and 47.86 for rotations per second, which is 
frequency . 
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cross-section 
turbine assent) I y 


Fig.1 Dimensions of Turl 





Turbine Housing 



Fig. 2 A Control Voume and a Schematic 
Diagram of a Turbine-Brush Unit 
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OUTLET 



(b)Combined Vector Diagram at Inlet, Exit 
of a Blade of a Turbine 


Fig. 5 Inlet and Exit Vector Diagram and a Combined Vector Diagram 


r2 dm = (density)(2*3.14rdr)(t) 



r2 


\ 


(b) Friction Force due to Centrifugal 
Force 


Fig. 6 Mass Moment of Inertia, Centrifugal Force, Friction Force 
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